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We study the possible heavy molecular states composed of a pair of charm mesons in the H and S 
doublets. Since the P-wave charm-strange mesons Dso(2317) and Dsi(2460) are extremely narrow, 
the future experimental observation of the possible heavy molecular states composed of Ds/D^ and 
-Dso(2317)/-Dsi(2460) may be feasible if they really exist. Especially the possible J^ = 1~~ states 
may be searched for via the initial state radiation technique. 
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I. INTRODUCTION 



The new family of the charmonium or charnionium-like states include X(3872), F(3940), y(4260), Z(3930), 
'^'. X(3940), r(4325), y(4360), r(4660), Z+(4430), Z+(4050), Z+(4250) and y(4140) etc [H-EI- Many states sit on the 
Qh' the threshold of two charmed mesons, which inspired some of them (especially those charged ones) to be candidates 
JL ,, of heavy moleculucs [12n20i | . 

rt\ • In the heavy quark limit, the S-wave and P-wave heavy mesons can be categorized into three doublets: H = (0^ , 1 ), 

^ ; S = (0+, 1+), T = (l+,2+). We collect their masses from PDG in Table HI The bottom mesons in the 5" doublet 
are still missing experimentally. Thus, we will adopt the theoretical predictions of the bottom meson masses in the S 
doublet when we study the heavy flavor molecular system composed of the bottom and anti-bottom mesons. 
'^ . In the framework of the meson exchange model, we have investigated the possible loosely bound molecular states 

T^ ' composed of a pair of heavy mesons in Refs. |16l - [2]| . In this work, we will investigate the possible heavy molecular 
0^ ■ system constructed by the charmed and anti-charmed mesons, where one meson is in the H doublet and the other 
0^ ' one is in the S doublet. In the following, we denote the heavy flavor molecular system as the HS + h.c. system for 
^^ I the convenience. 

l/~j ■ The HS + h.c. system can be categorized into four subsystems: [PPq], [PPi], [P*Po] and [P*P{]. They correspond 

^^ [ to different quantum number combinations 0~ -I- 0+, 0~ -I- 1"*", 1~ -I- 0"^ and 1~ -I- 1+, respectively. Since charmed 
^^ . mesons belong to the fundamental representation of flavor SU{3), the system constructed by the charmed meson and 
anti-charmed meson forms an octet and a singlet: 3(8)3 = 8®las illustrated in Table HIl The parameter c = =Fl 
in the flavor wave functions corresponds to the charge parity C ~ ±1 respectively for the neutral systems as pointed 
^. out in Refs. [iMl. 

rS I This paper is organized as follows. We review the formalism in Section |lT] and present the results in Section IIIII 

j^ ' The last section is a short summary. 

II. THEORETICAL FRAMEWORK 

A. The potential model 

The potential model is an effective approach to study the two-body bound state problem. For the HS + h.c. system, 
the scattering between the charmed and anti-charmed meson occurs via exchanging the light pseudoscalar, scalar and 
vector mesons, which play the role of providing long-distant, intermediate-distance and short-distance forces. At the 
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hadron level, there exist two types of diagrams in the scattering of the charmed and anti-charmed mesons, i.e., the 
cross and direct diagrams as shown in Table Hill The exchanged light mesons relevant to the four subsystems [i^i^g*]' 
[PP{], [P*Po] and [P*P{] are also presented in Table lllll When writing out the scattering amplitude, the monopole 
form factor is introduced at every interaction vertex to compensate the off-shell effect of the exchanged light meson 



-^(9) 



A2 



m 



A2 



(1) 



where the phenomenological cutoff parameter A is about 1 GeV. q and m denote the four-momentum and the mass 
of the exchanged meson. 







H doublet 








charmed-up(down) meson 


charmed-strange meson 


J^ 


charged 


mass (MeV) neutral mass (MeV) 


charged mass (MeV) 


0~ 

1" 




1869.3 D° 1864.5 
2010.0 D*° 2006.7 


Dt 1968.2 
D:± 2112.0 






bottom-up (down) meson 


bottom-strange meson 


J^ 


charged 


mass (MeV) neutral mass (MeV) 


neutral mass (MeV) 


0" 

1" 


B± 
B*± 


5279.1 B° 5279.5 
5325.1 B*° 5325.1 


B° 5365.1 
Bf 5412.0 



S doublet 







charmed-up(down) meson 


charmed-strange meson 


J^ 


charged 


mass (MeV) neutral mass (MeV) 


charged mass (MeV) 


0+ 

1+ 




2403 Dg 2352 
Df 2427 


D% 2317.3 
D'^ 2458 






bottom-up (down) meson 


bottom-strange meson 


J^ 


charged 


mass (MeV) neutral mass (MeV) 


neutral mass (MeV) 


0+ 

1+ 




B? 


Bso 
Bs\ 



TABLE I: The masses of the heavy mesons in the H and S doublets [22 



According to the effective Lagrangian describing the interaction between the light and charmed mesons, one can 
write down the scattering amplitude between the charmed and anti-charmed mesons. Such a system can be described 
as 



Ai ,A2 



(2) 



where J and Ji [i = 1,2) denote the angular momentum of the system and the componentents. The scattering 
amplitude iAi{J, Jz) is related to the interaction potential in the momentum space in terms of the Breit approximation 



V(q) = - 



1 



n, 2m,^f2m/ 



=X(J,Jz) 



(3) 



Here, mi and rrij denote the masses of the initial and final statrs respectively. The potential in the coordinate space 
V(r) is obtained after Fourier transformation. 




*]"® 




Jsl 



3^3=801 



(a) [PPo*] system 


States 


Wave function 


[PPS]+ 


^^[D+D?, + cD'Dt] 


[PPof 


i[(D°5S+cO«D«) 




-{D+D-+CD-D+)] 


[ppsr 


^[D-D?, + c&'D^] 


[PPS]t 


^[D+5« + cD"D+] 


[PPofs 


^^[DtD-+cD-Dto] 


[PPS]s 


^[D-D?,+cD"D-,] 


[ppsn 


j=[D+D:, + cD-D+] 


[PPofs 


i[(D"5!J+cO°Dg) 




+iD+D-+cD~D+)] 


[PPofk 


^[DtD;, + cD:Dto] 



(b) [PPi] system 


States 


Wave function 


[PPi] + 


^[D+DT + cD'>D'+] 


[ppir 


^[{D"Df + cD°Df) 




~{D+D^+cD-D+)] 


[ppi]- 


-^[D-DT + c&'D'f] 


[ppi]t 


^^[DtD[" + c&'D'+] 


[ppiTs 


^^[DtD[-+cD-D'+] 


lppi]J 


^^[D-DT + cD'D'-] 


[ppin 


^[D+D'-,+cD7D'+] 


[ppi]'^ 


\[{D"Df + cD°D'i') 




+(d+d;-+cd-d'+)] 


[PPifsl 


A^[DtD'r,+cD-D'+] 



(c) [P*PS] system 



(d) [P*Pi'] system 



States 



[P'PS]+ 
[P'PSV 

[p'psr 
[p*Po]t 
[p'Pot 

[P*Po*]9 
[P*Po1'^ 

[P*PSti 



Wave function 



States 



^{D'+ D^ + cD- D+)] 

^^[D*- 01 + 00*^0^] [P*Pi]- 

^[D:+D°o+cD"'Dto] [P*Pi]t 

^[D:+D^+cD'-Dto] [P'Pifs 

^^[D:-D°+cD"'D:,] [P'Pi]- 

^[D'+d;, + cd:-d+] [P*A']9 

^[D:+D-^+cDrDto] lP''Pi]°i 



Wave function 



i[(P>*Op)f + ££»*» Df) 
-{D'+D-+cD*-D+)] 
^[D'-D'° + cD*''D'-] 
^[Dl+D? + cD*°D'+] 
^[D:+D'{-+cD*-D'+] 
^^[D:-DT + cD*^D'-,\ 
^[D*+ D'-, + cDT D'+] 

+{0*+ D'^- + cD- D'+)] 
^[Dl+ D'- + cDT D'A\ 



TABLE II: The flavor wave functions of the [PPq], [PPi], [P* Po] and [P*Pi] systems. We use [■ ■ ■]l to name the different 
states, where the superscript a denotes the charge of the state and subscript b is introduced to distinguish the different states 
in a subsystem. 



B. The effective Lagrangian 



The effective Lagrangian describing the interaction between the hght and heavy flavor mesons is constructed with 
the help of the chiral symmetry and heavy quark symmetry [23| 

C = igTr[Hb-i^At^^5Ha] + ig'Tr[Sblt,Al^l^Sa] + {ihTrlSb-t^A^.l^Ha] + h.c.} 
+iPTr[HbV^{V^ - Pf,)baHa] + iXTr[Hba''''F^4p)baHa] 
+i/3'Tr[Sbv''{V^ - p^)baSa] + iX'Tr[Sba^''F^,{p)i,aSa] 
+{iCTr[Hb-f^{V^ ~ P^)6a5j + ipTr\HbO>'''F^,{p)haSa\ + /i-c.} 

^,Tr\UaGH,\+ g';rr\Sa<jSa\ + \^TT\Sa{d^,a)^^Ha\ + h.c\. 



~9a 



The H and S fields, which correspond to the (0,1) and (0+, 1+) doublets are defined respectively 

1 



ff„ 



■K^T'^+^^aTs], 



Hb = -/'H^i" = [p:ir+'^pU]^, 



Sa = 



^-^[Kl,l5 - Poah 



S, = 7°^^° = [Pi'L7^75-PoJ] 



*til + ^ 



(4) 



(5) 



where v = (1,0,0,0). Here, the annihilation operations P, P* Pq, P{^ are of dimension 3/2 and satisfy the normal- 
ization relations 



(0|P|gq-(0-)) = ,/Mh, (0|P;|Qg(0-)) = c^^Mh, 
(0|Po*|gg-(0+)) = ^/Ws, {Q\P[^\Qq{l+)) = e^./Ws. 

In Eq. ((4]), the expansion of the axial vector gives 



(6) 



with ^ = exp(iP//^) and /^ = 132 McV. p^^ = ig^.W^^/V^ and F^,^{p) = df^pi, - d^p^ + [p^, p^]. The octet 
pseudoscalar and nonet vector meson matrices read as 



(7) 





7:+ K+ ^ 






p+ K*+ ^ 


n^ 


-S + ^ ^" 


, V = 


p^ 




[ ^- 


^° -fj 




I ^*- 


K*^ J 



After expanding Eq. Q , the effective Lagrangian of the pseudoscalar mesons with heavy flavor mesons reads 

/:i,.OT = ^'^{PbP:l-p:^Pl)^^na, 

Ct>t,iy = ^{2^v;)(Pl,Pl-PbP;l)^^^ba. 
Cv'V'.v = i^{P[b ■ p:^ + p: ■ Pii)i2iv) ■ dFba ■ 



(8) 

(9) 

(10) 

(11) 
(12) 







Direct diagram 




Crossed diagram 




A 


Pi ->- 


^ 


J 


Pi ->- 


P3 -^ 








^^t 


Exchanged 






Q T Exchanged 










mesons 








mesons 




B- 




1 1 n 


B- 




ii T 




po -^ 


PI ->- 


P2 ->- 


Pi ->- 


Subsystems 


Pseudoscalar Vector Scalar 


Pseudoscalar Vector Scalar 






(P) 


(V) (§) 




(P) 


(V) (§) 


[PPol 






/ / 




/ 




[PPi] 






/ / 






/ / 


[P'PS] 






/ / 






/ / 


[P'Pi] 




/ 


/ / 




/ 


/ / 



TABLE III: The direct and crossed scattering diagrams with the exchanged hght mesons for the HS + h.c. system 



The effective Lagrangian of the vector mesons with heavy flavor mesons reads 



C 



vv 



;^9Vp* p*U 



,^9\ 



-i^P: ■ P:\2iv ■ Yt,a) - ^4:^P;^Pr^(5MV. - d,Y^),a, 



Cv'„v--v = -i—i^PobPoli2iv -Wba), 



C 



vxv 



,C5 



^3V(0n.,a\(p*Pp*t 



2^-^{p;^p;l + PSbPr^)Y,ba 2^{2^v-){p;^p;l p^bP*J^){d^Y p a^v„) 



V2'^' 
9v^ 

/2 



V2 



ha 1 



r^j^iv = ^^Plfc•P{i(2^«.Vfc,) + ^4^PlVP;L(a^V'^-a''V'^),,, 



,^5v, 



/:p;PV = ^2^(nPr - P[tPl)N,ba + ^2^(2^«")(nP;f + P^f P2)(a„V, - dpN^U. 



/:^jp.v = ^(2zt,^)£.,^,(prpf; - p'^bPr^)^^ + 2^£.„^^(prpr + Pi'^pr^)(9"v^ - 9^v"),„, 

£^Piv = ^2^(nP;;:t - PiPit)V,,, - ^2^{2^vP){PbP[Z^ + P;?P+)(a„V, - 9,V.),., 
£p;pjv = ^2^(2^^;^)e,„^^(Po;p{;:t + P(^Po*i)(5"V^),,. 
The effective Lagrangian of the a meson interacting with heavy flavor mesons reads 

CvVa = -2go-0-P;,P/, 

Cv'V. = 2g^aP:-P*\ 
C-v-V'^a = 2g'„aP^^P*l, 



(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
(22) 

(23) 
(24) 
(25) 
(26) 
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Ct^t^',. = -^^{Pi,,P^ - PbP[l)d''a, (27) 






Cv'V',. = -'^-^{PmP:'^^P:^P;l)dxa, (28) 

Cv'v-a = ^^^v'^eo.,^APTaP^^ -P:''P[':)d^'y- (29) 

The values of the parameters g'^'\ (3'^'\ A''', ft., zu, C, g^ and /lo- arc discussed in Refs. |23|.|24|. 

C. The effective potential of the HS + h.c. system 

We list the effective potentials of the states in the HS + h.c. system in Tabic ITVl Here V^ , Ql, W^ and O^ are 
the sub-potentials corresponding to the subsystems [PP^], [PP[], [P*P^] and [P*P[] respectively. The parameters 
in front of the sub-potentials in Tabic IIVI are from the coefficients in the octet pscudoscalar and nonet vector meson 
matrices. 

For the [-P-Pq ] subsystem, the sub-potentials are 



yD^rect^^■^ = 


~-liP'glY{K,qo,m^,r), 


ymrect^^-^ = 


--9<jgaYiA,qQ,m„,r), 


l^P^™-(r) = 




r the [PP[] subsystem, the sub-potentials are 




Q^''-'==*(r) = 


-^l3^'gvY{A,qo,mv,r), 


QD^rect^^^^ = 


-9<7gaY{A,qo,ma,r), 



(30) 

(31) 

(32) 



(33) 
(34) 



qC™..(^) = -^i2e~8Cwq'o + 8^'q'o^)Y{A,qim^,r) 



-f(8.^-^ 
12 m 



{8fi^-Z^)Z{A,q'o,my,r), (35) 



The sub-potentials relevant to the [P*Pq] subsystem are 



QCross^^^ = -ilz(A,g^,™.,r). (36) 



Wmrect^^^ = -^(3l3'glY{A,qo,my,r), (37) 

W^^^^^\r) = -g,g',Y{A,qo,m,,r), (38) 
W^-^'^r) = ig2,(2C2-8Ctng^ + 8n72g^')r(A,q^,mv,r) 

-^gvi8^J-^~%)Z{A,q'o,m^,r), (39) 

M^^c.— (r) = ^lziA,qim^,r). (40) 



States 


Effective potential 


[pp^.]±.0 


[ 


_iymr.ct^^-^ + lV^^^---\r) + Vf'^'^'ir)] + c[-iyf --(r) + ^V^^^'^'ir)] 


[ppl]±,0 


h^Q? 


"■^^*(r-) + iO^''-^^*(r) + Q?"-^^*(r)] + c[-iQ^--(r) + iQ5'-°^=(r) + 0?'-°^^(r)] 


[p*p^.]±,0 


[- 1 Vl/D--*(r) + iM/i'''--'(r) + W^''-'"'\a)] + c[~-^W;^^''"'ir) + lwS''°''{r) + W^'''"'ir)] 


[p.p^/]±,0 






[PPo]t'", [PPofs 




c[-il/^c— M] 


[PP{]±'0, [PPi]° 







jp*p^.j±.o^ [P'Po-]« 







[P*Pl']±'0, [P'Pi'l? 




[-10f-*W] + chiO,f— (r)] 


[pporsi 




[n'''""''W] + c[|Kf'°=^(0] 


[PPifsl 




[Qf-*(r)]+c[Q^'-— (r)] 


[P*Pofsl 




[VKf'"=*(r)]+c[VK^^'-°^=(r)] 


[P'Pifsl 




[Of ==*(r) + f Of ^'=^*(r)] + c[Of °-(r) + |0^'-°^^(r)] 


lPPS]°s 




[1^/?— '(^) + i^f "■-^(r) + l//'''--*(r)] + c[|Ff '■°-(r) + ^^^^"^^(r)] 


[PPifs 


ilQ? 


'"='(0 + iO^''-'=='(r) + 0?''-^='(r)] + c[|Q^'-°-(r) + iQS'-°^=(r) + O^'-^(r)] 


[P'Pofs 


[3y^,mrect^^■^ + ^W^'^'^^^r) + W^'^^^\r)] + c[|Wp^^°=^(r) + iM^^'-°-(r) + W^^''''{r)] 


[P*pi]<i 




[|0^''-==*(r) + i05"""'(r) + fOf "■•=^*(r) + iOf '=^*(r) + 0?'^'=^*(r)] 



TABLE IV: The total effective potentials of the states in the subsystems [PPq], [PPi], [P*Po] and [P*Pi'] 



For the [P*P[\ subsystem, the relevant potentials are 

1 



25^/3/?'C[J]y (A, qo, mv, r) - 2gl\\' B[J]Z{K, qo,mv, r), 
of ^^*(r) = -9,5;0[J]y(A,go,mv,r), 






-A\J\Z{k,qa,my,r), 



o: 



Cross 



1 



(r) = -g^y{2C^^8Cwq'o + 8m^q'„')C[J]Y{A,q'o,mv,r) 

1 2(2 
4 m^ 

0-°-(r) = -^q',^C[J]YiA,qimr,r), 



(41) 
(42) 
(43) 

(44) 
(45) 
(46) 



with the coefficients A[J], B[J] and C[J] 





J=0 


J=l 


J=2 


A[J] 
B[J] 
C[J] 


2/3 

4/3 

1 


1/3 

2/3 

1 


-1/3 

-2/3 
1 



(47) 



In the above expressions, the functions Y{A, k, m, r) and Z{A, k, m, r) are defined as 



y(A, K,m,r) 



_A_ iq.r 

(27r)3 

if IkI < m 



-m? \ A^ — q2 



(48) 



_ otherwise, -^^ (cos(CIr) - e-?^-^) + ^ ((2 + §) e 



<2r 



and 



Z(A, K, TO, r) 



1 a 
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r — -y(A, K, 771, r) 
or 



(49) 



(70 = U and ttie mass differences q'^ as 



where Ci = ^/m? — k^, C( = Vk^ — tti^, and C2 = x/A^ — k^. We take (70 = and the maoo umcicii^co yg txo 
(mpo - TOp), (top^' - TOp), (top„ - TOp.) and (mp^' - mp.) corresponding to the [PPq], [PPi\, [P*Pa] ^nd [P* P[] 
subsystems respectively. 



III. NUMERICAL RESULTS 

In the following, we present the numerical results for the different systems. Throughout this work, jjl denotes the 
reduced mass of the corresponding system in all figures of the numerical results. 



A. The [PPS] system 



[PPS] 



The total potentials of the [PPq] '^ system with c = ±1 and the typical values of A arc presented in Fig. [TJ The 
comparison of the total potential with the partial potentials indicates that the tt exchange plays the dominant role 
in the total effective potential of the [PPq ]^'" system with c — ±1. Here we only illustrate the potential with A = 1 
GeV in Fig. [1] (b) and (d). 

Since the one tt exchange potential is dominant in the total potential of [PPJ]^'", it's enough to consider the one 
TT exchange force only when we study the bound state solution of the [PP^]^'^ system. The one tt exchange potential 
is proportional to the coupling constant h. Thus, we try to solve the schrodinger equation with the obtained one pion 
exchange potential under some typical values of h. 

For the [PPg*]^'° system with c ~ +1, we can find the bound state solutions numerically. As shown in Fig. [21 
the binding energy Eq decreases when A becomes larger. One may further check the corresponding wave function to 
make sure whether the bound state solution is reasonable. An example is shown in Fig. [2] (a) with h = —0.56. As 
A increases, the wave function oscillates at large distance. In other words, smaller A values such as A < 0.75 GeV 
lead to the bound state wave function with less oscillating behavior. On the other hand, the reasonable range of A is 
sometimes assumed to be around 1 GeV. 

Enhancing the coupling constant helps to form the bound state of [PPq*]='='° system with c = +1. For example, 
there exist a very reasonable bound state solution with 1.3/i corresponding to a shallow binding energy, a reasonable 
A and wave function with good behavior as shown in Fig. [3] (b). If the coupling constant becomes even larger, the 
binding energy of the [PP^]^'^ system with c = +1 becomes deeper. As shown in Fig [5] (c) and (d), there only exist 
deeply bound states for the [PP^]^'^ system (c — +1) with 2h and 3/i. 



For the [PPq]^'° system with c = —1, we can not find the bound state solution with Ih. With 2h and A = 2.6 ~ 3 
GeV, the wave function of the [PPq]^'^ system with c = — 1 becomes non-oscillating. When further increasing the 
coupling constant to 3h, there exists a [PPq]^'^ bound state with c = — 1 and A near 1 GcV. The numerical results 
of the [PPq]^'" system with c = — 1 are presented in Fig. [3] 

2. [PPo*]?'", [^^oli 

The potential of [PPq*]^'", [PPq]^ is shown in Fig. 01 where only the rj meson exchange contributes. The effective 
potential is repulsive and attractive for c = +1 and c = — 1 respectively. There does not exist any bound states with 
c = +l. 

For the [PPg]^'", [PPo]| system with c = —1, we can not find the bound state solution with h = —0.56. We 
illustrate the numerical results with 2h and 3h in Fig. [5l There may exist a [PPg*]^'°, [PPg*]? bound state with 
c = — 1 if the coupling constant is around several h. 

3. [PPSti 

In Fig. [HI we plot the variation of the effective potential of the [PPo*]si system with r and different coupling 
constants and the cutoff parameter. The rj and (j) meson exchange contributes to the total effective potential. The 
cutoff A should be larger than the mass of (f) meson. The total effective potentials with (c = +l,/3/3' < 0) and 
(c = —1,(3/3' < 0) are attractive while the effective potentials with (c = +1,(3/3' > 0) and (c = —1,(3/3' > 0) are 
repulsive. The comparison of the total potential with the ?; and meson exchange potential is given in Fig. [5] (e) and 
(f) with A = 1.5 GeV. 

For the [PPoJ^i system with c — —1, there docs not exist a bound state not only for the case (3/3' > but also for 
(3(3' < 0. For the [PPq ]si system with c = +1 and (3(3' > 0, a reasonable solution of bound state appears with 3h. For 
the [PPq ]si system with c = +1 and /3/3' < 0, there does not exist a bound state with Ih. The bound state solution 
appears with 2h and 3h. The detailed numerical results of the [PPo*]si system with c = +1 arc shown in Fig. [T] 

4- [PPofs 

The dependence of the total potential of the [PPo*]8 system on A is presented in Fig. [51 The total effective potential 
oscillates with r. The n, rj, a, p and a; exchange contributes to the total exchange potential. With A = 1 GeV, the 
corresponding tt, 77, a, p and a; exchange potentials for the [PPo*]g systems with c = ±1 are given in|Hl(c) and (d). The 
TT exchange potential is dominant. Therefore, we only consider the single pion exchange potential when we explore 
whether there exists the bound state solution of the [PPq ]§ system. 

For the [PPqJs system with c = +1, there exists a bound state solution. The variation of the binding energy of 
[PPpJg with c — +1 is shown in Fig. [51 with different values of the coupling constant. The wave functions is not 
reasonable when taking \h\ = 0.56, 0.65, 0.7, which corresponds to the shallow binding energy of the [PPqIs system 
with c = +1. Increasing the coupling constant to 2h leads to a bound state solution with a reasonable wave function. 

The numerical result presented in Fig. [TOl shows that there exists a [PPo*]8 bound state with c = — 1. 

B. The [PPi'] system 
1. [PPi]^'"' 

The a, p and lu meson exchange contributes to the total effective potential of the [PP{]='='*^ system. Since the p and 
Lo exchange potentials cancel each other, the a meson exchange potential is dominant. We only need to consider the a 
exchange potential in order to explore whether there exists a [PP;[]^'° bound state. In Fig. [TTl the effective potential 
of the [PP!^'^''^ system is presented. The a exchange potential of the [PP^']^'*' system with c = +1 is repulsive. There 
does not exist a [PP^^^'^ bound state with c = +1. 

For the [PP'^\^''^ system with c = — 1, the a exchange potential is attractive. Meanwhile, the cross diagram plays 
the dominant role. However, in the cases Qag'^ > and §09'^ < 0, we can not find a bound state solution of the 
[PP[\^'^ system with c = — 1 and the coupling constants listed in the caption of Fig. [TT] when scanning the range of 
cutoff A. 



10 

When taking 2/icr, we present the variation of the binding energy of the [PPl]^'^ system with c — ~l and the cutoff 
A in Fig. [T2]in the cases ga9a > and ga9a < 0. However, the corresponding cutoff A is far larger than 1 GeV. Thus, 
we tend to conclude that there does not exist a [PP[\^'^ bound state with c = —1. 

2. [PPi]t-°, [PPif, 

There does not exist the [PP^]^'", [PP[\^ bound state system since none of the meson exchange forces is allowed 
as indicated in Table HVl 

3. [PPiTi 

Here the effective potential of the [PPi]^i arises from the (j) meson exchange. There exist eight combinations of the 
parameters c, (3/3' and (m. With different parameters, the dependence of the effective potential of the [PPi]si system 
with c = ±1 on r is shown in Figs. [13lfT4l 

When solving the Schrodinger equation, the cutoff A is expected to be larger than the mass of the (f> meson. We 
only find bound state solutions of [PPiJ^i with c — ~l, f3(3' < and d^w < 0, which is shown in Fig. [15] 

4- [PPi]t 

The a, p, to exchange contributes to the total effective potential of [PPilg system. In Fig. [161 the variation of the 
total effective potential of [PPi'Jg to r and different typical A values is given. The comparison of the total potential 
with the a, p and uj exchange potential is presented in Fig. [THJ (c) and (d) with A = 1 GeV. 

Since the total potential is repulsive for the case of c = +1, we only explore wether there exists the bound state 
[PP{]g system with c = —1. The combinations of coupling constants go-g^i I^P' ^-nd C'^ a-re eight. However, we find 
the bound state solution of [PP{]g system with c = — 1 with six parameter combinations just listed in the caption of 
Fig. [HI 

C. The [P*Po] system 

1. [P*Po*]±'° 



The total effective potential of the [P*Pg*]*'° system arises from the a, p and to meson exchange. Since the p 
and id exchange potentials cancel each other, the a exchange potential is dominant in the total effective potential 
of [P*Pg*]^'''. We only consider the ui exchange potential contribution to explore whether there exists a [P*Pg*]^'° 
bound state . In Fig. [TBI the effective potential of the [P*Pg*]^'° system is presented. 

The a exchange potential of [PP[\^'^ system with c = — 1 is repulsive. Thus there does not exist a [PP{]^'^ bound 
state with c = — 1. 

For the [P*Pg*]^'° system with c = +1, the a exchange potential is attractive. Meanwhile, the cross diagram plays 
the dominant role to the a exchange potential. However, in the gag'a > and g^g'^ < two cases, we can not find 
bound state solutions of the [P*Pq]^'° system with c = — 1 and the coupling constants listed in the caption of Fig. 
[iHJwhen scanning the range of the cutoff A. 

When taking 2ha, we obtain the binding energy of the [P*Pg*]='='° system with c = +1 dependent on the cutoff A 
in Fig. [T9|only in the gag'a- < case. However, the corresponding cutoff A is far larger than usual 1 GeV. Thus, the 
[P*Po*]^'° system with c = +1 seems impossible to form the bound state. 

2. [p*Po*]^^ [p*Po]i 

There does not exist the [P*Pq]^'*', [P*Po*]° bound state since no suitable meson exchange force is allowed as 
indicated in Table ITVl 
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3. [P*Po*]°i 

Here, the effective potential of the [-P*i^o*]si system resuhs from the (j) meson exchange. There exist eight combi- 
nations of the parameters c, /?/?' and C^vj. Under different parameter space, the dependence of the effective potential 
of [P*Po1°i with c = ±1 on r is listed in Figs. HDHH] 

When solving the Schrddinger equation, the cutoff A should be larger than the mass of the meson. There exists 
a bound state solution of the \P*PQ]^gi system only with c = +1, /3/3' < and C,tu < 0, which is shown in Fig. [521 

4. [P'Pofs 

As indicated in Table ITVl the behavior of the p meson exchange potential is similar to that of the uj meson exchange 
potential. When exploring the [^*-Po ]§ system, we need to consider the a exchange potential together with both the 
p and Lo meson exchanges. The relevant total and partial potentials are shown in Fig. 1231 

Since the total potential of the [i-'*Pg*]g system with c = — 1 is repulsive, there does not exist a [P*Po]8 bound state 
system with c = —1. In the following, we mainly focus on the [i^*-Po]8 system with c = +1, where the total effective 
potential is attractive. As shown in Fig. [511 there exists the bound state solution of the [P*Po Jg system with c = +1. 

D. The [P*Pi'] system 

In Rcfs. [ll,[il, [H], the [P*P{]'^'^ and [P*Pi]s systems were studied. In this work we will discuss the rest two 
cases of the [P*Pi] system, i.e., [P*Pi]t'°/[P*Pi]°s and [P*PXi- 

1. [P'Pi]f'°,[P'Pi]'i 

The effective potential of the [P*Pi]f''^, [P*Pi]^ system arises from the t] meson exchange only, which is plotted in 
Fig. 1251 under several typical combinations of parameters. 

When solving the the Schrodinger equation, we obtain the bound state solutions of the [D*D'i]f'^, [D*D'i]^ system, 
which are listed in Fig. [251 The bound state solution of the [D*D[]f''^, [D*D[]'-g system appears when (c = +1, J = 
0, gg' > 0), (c = +1, J = 0, gg' < 0), (c = +1, J = 1, gg' > 0), (c = +1, J = 1, gg' < 0), (c = +1, J = 2, gg' > 0) 
and (c = +1, J = 0, gg' < 0). The above numerical results are obtained with 3/i while we can not find the bound 
state solution for Ih and 2h. 

2. [P''P^°si 

The effective potential of the [P*P{]2i system results from both the i] and 4> meson exchanges as shown in Figs. [27l 
andm 

There exist 96 combinations of the parameters c, /3/3', gg' ^ XX' and (^w. When solving the Schrodinger equation, 
the cutoff A should be larger than the mass of the (j) meson. We only find bound state solutions of the [P*Pi]'^i 
system under 20 combinations of c, 13/3', gg', XX' and (zj, which are shown in Figs. [^ 1501 and [5l1 for J = 0, 1, 2 cases 
respectively. 

IV. CONCLUSION 

We have investigated the possible molecular states composed of a pair of heavy mesons in the H and S doublet. 
They are expected to be loosely bound states mainly via the long range pscudoscalar scalar meson exchange. We are 
interested in these shallow S-wavc states with J^ = 0~, 1~,2~, especially those neutral ones with cither C = + or 
C = — . Those states with J^*^ = 0^^, 1 ^ are exotic. They may exist with very reasonable coupling constants. For 
example, the J^'^ = state is shown in Fig. [2l 

The non-strange P-wave (0+,l+) heavy mesons are very broad with a width around several hundred MeV [22 |. 
Instead of forming a stable molecular state, the system containing one non-strange P-wave heavy meson may decay 
rapidly. Experimental identification of such a molecular state may be difficult. The attractive interaction between 
the meson pair may lead to a possible threshold enhancement in the production cross section. 
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In contrast, the P-wave charm-strange mesons Dso{2317) and I?si(2460) lie below the DK and D*K threshold 
respectively. They are extremely narrow because their strong decays violate the isospin symmetry. The future 
experimental observation of the possible heavy molecular states composed of Ds/D* and Dso(2317)/-Dsi(2460) may 
be feasible if they really exist. 

There may exist (1) two 0^ states around 4.25 GeV with different C-parity; (2) two 0^ states around 4.5 GcV 
with different C-parity; (3) two 2" states around 4.5 GeV with different C-parity; (4) four 1^ states around 4.35 GcV 
with different C-parity; (5) two 1~ states around 4.5 GeV with different C-parity. The three neutral J^*-^ = 1 
states may be searched for via the initial state radiation technique. Wc notice that they are are rather close to those 
1~~ states observed by Belle around this mass region. The other states might be produced from B or Bg decays if 
kinematically allowed or at other possible facilities such as RHIC, Tevatron and LHCb. 

The dominant decay modes of the above states are the open-charm modes Dg Ds if angular momentum con- 
servation, parity and C-parity symmetry allow. The other characteristic decay modes are the hidden-charm modes 
containing one J/tp, or ijc or Xcj etc. One may easily exhaust the possible final states according to C/P parity and 
angular momentum conservation and kinematical considerations. For those non-exotic states, they may be signifi- 
cantly narrower than the conventional charmonium around the same mass region because of their molecular nature. 
The manifestly exotic states may be searched for through their quantum number. 
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The total potential function 
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FIG. 1: (a) The variation of tiie total potential of the [DDo] system (c = +1) with r and A; (b) the exchange potentials 
of the TV, T], a, p, uj mesons of the [DDo] system (c = +1) with A = 1 GeV; (c) the variation of the total potential of the 
[DDo] system (c — —1) with r and A; (d) the exchange potentials of the vr, 77, a, p, ui mesons of the [DDo] system 
(c = —1) with A = 1 GeV. The above potentials are obtained when g^g',^ > and /?/?' > 0. Here, \h\ — 0.56, \ga\ = 0.76, 
\g'^\ = 0.76, |/J| = 0.909 and |/3'| = 0.533. 
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FIG. 2: The dependence of the binding energy of the [DDq] system (c = -1-1) on A and the wave function of this system. 
Here, we only consider the single n exchange potential, (a), (b), (c) and (d) correspond to the case of l/i, 1.3/i, 2h and 3/i with 
\h\ = 0.56 respectively. 
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FIG. 3: The dependence of the binding energy of [DDo] ' system (c = — 1) on A and the wave function of this system. Here, 
we only consider the single n exchange potential, (a) and (b) correspond to the case of 2h and 3/i with \h\ = 0.56 respectively. 
When taking Ih, we can not find the bound state solution. 
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FIG. 4: (a) The variation of the total potential of the [DDoJ^ ''\ [DDo]i system (c = -1-1) with r and A; (b) The variation of 
the total potential of the [DDo]s ' , [DDo]g system (c — —I) with r and A. Here, we take \h\ — 0.56. 
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FIG. 9: The dependence of the binding energy of tlie [DDoJs system (c — +1) on A and the wave function of this system. 
Here, we only consider the contribution of the tt exchange potentials. We take different values of h, in diagrams (a)-(e). The 
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FIG. 11: (a) The variation of the total potential of the [DDi] ' system (c = +1) with r and A; (b) the variation of the total 
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FIG. 13: 
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(a), (b), (c) and (d) are the potentials ol the [PP[]°i with (c = +1, ^/3' > 0), (c = +1, ^/3' < 0), (c = -1, ^/3' > 0) 
-1, PP' < 0), respectively. Here, we take C^ > with |/3| = 0.909, |/?'| = 0.533, |C| = 0.727 and \vj\ = 0.364. 
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FIG. 14: (a), (b), (c) and (d) are the potentials ol the [PP[]°i with (c = +1, PP' > 0), (c = +1, PP' < 0), (c = -1, PP' > 0) 
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FIG. 16: (a) and (b) are the total potentials of the [PPi]s system with c — +1 and c = — 1 respectively when gag'cr > 0, /3/3' > 
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FIG. 17: The bound state solution of the [PPi']g with c — —1. (a), (b), (c), (d), (e), and (f) correspond to the results with 
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(g'rg'cr < 0, /3/3' >0, Cu7> 0), we can not find the bound state sohition. Here, \g„\ = 0.76, \g'„\ = 0.76, \h„\ = 0.323, |/3| = 0.909 
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FIG. 18: (a) The variation of the total potential of the [D*I?o] system (c = -1-1) with r and A; (b) the variation of the total 
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FIG. 20: (a), (b), (c) and (d) are the potentials of the [D*-D5]si with (c = +1, /5/3' > 0), (c = +1, ^/3' < 0), (c = -1, /3^' > 0) 
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FIG. 21: (a), (b), (c) and (d) are the potentials of the [D* D^fs^ system with (c = +1, /3/3' > 0), (c = +1, pp' < 0), 
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FIG. 22: The dependence of the binding energy of [D*A*]°i on A with c = +1, /3^' < and (cc < with j/Jj = 0.909, 
\I3'\ = 0.533, Id = 0.727 and \zj\ = 0.364. 
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FIG. 24: The bound state solution of the [D*DQ]g with c = +1. (a), (b), (c), (d), (e), and (f) correspond to the results with 
{g.g'^ > 0, 13 p' < 0, Ct^ > 0), {g.g'^ > 0, /J/3' < 0, C^ < 0), {g.g', > 0, /3/3' > 0, C^ < 0), {g.g', < 0, /3/?' < 0, C^^ > 0), 
(ffcTpCT < 0, /3/3' > 0, Cuj < 0) and {g^ga < 0, /3/3' < 0, C^^ < 0) respectively. When taking {g^ga > 0, fi/S' > 0, (uj > 0) or 
(,g<T5a < 0, /3/?' > 0, C^ > 0), we can not find the bound state solution. Here, \ga\ = 0.76, \g'^\ = 0.76, \h^\ = 0.323, |/J| = 0.909 
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FIG. 25: The variation of the potential of the [D* D'l]^ ' , [D* D'i]g system, where we only list several typical cases with 
parameters (c = +1, J = 0, gg' > 0), (c = -1, J = 0, gg' > 0), (c = +1, J = 1, gg' > 0), (c = -1, J = 1, gg' > 0), 
(c = +1, J = 2, gg' > 0) and (c = —1, J = 0, gg' > 0), which correspond to diagrams (a), (b), (c), (d), (e) and (f) respectively. 
The above results are obtained with Ih. Here, \h\ = 0.56, \g\ = 0.75 and \g'\ = 0.25. 
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FIG. 26: The bound state solutions of the [D*Di]^'°, [D*D[fg system with (c = +1, J = 0, gg' > 0), (c = +1, J = 0, gg' < 0), 
(c = +1, J = 1, gg' > 0), (c = +1, J = 1, gg' < 0), (c = +1, J = 2, gg' > 0) and (c = +1, J = 0, gg' < 0), which correspond 
to diagrams (a), (b), (c), (d), (e) and (f) respectively. The above numerical results are obtained with 3h while we can not find 
the bound state solution for l/i and 2h. Here, \h\ = 0.56, \g\ = 0.75 and \g'\ = 0.25. 
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FIG. 27: The variation of the potential of the [D* D'i]^i system with r. (a), (b), (c), (d), (e) and (f) are for parameters 
(c = +1, J = 0, PP' > 0, AA' > 0, C^ > 0), (c = -1, J = 0, /3/3' > 0, AA' > 0, C^ > 0), (c = +1, J = 1, /3/3' > 
0, AA' > 0, Ccc7 > 0), (c = -1, J = 1, /3/3' > 0, AA' > 0, C^ > 0), (c = +1, J = 2, /3/3' > 0, AA' > 0, C^ > 0) and 
(c = -1, J = 2, /?/?' > 0, AA' > 0, C^ > 0), respectively. Here, \h\ = 0.56, \g\ = 0.75, \g'\ = 0.25, |A| = 0.4546, |A'| = 0.2667, 
|/3| = 0.909, \I3'\ = 0.533, |C| = 0.727 and |t;<7| = 0.364. 
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FIG. 28: The comparison of the rj and exchange potentials with the total effective potential for the [D*D'i]'^i system with 
A = 1.5 GeV. (a), (b), (c), (d), (e) and (f) are for parameters (c = +1, J = 0, /3/3' > 0, AA' > 0, gg' > 0, (zu > 0), 
(c = -1, J = 0, PP' > 0, AA' > 0, gg' > 0, C^ > 0), (c = +1, J = 1, pp' > 0, AA' > 0, gg' > 0, (uj > 0), (c = -1, J = 
1, PP' > 0, AA' > 0, gg' > 0, (zu > 0), (c = +1, J = 2, /J/?' > 0, AA' > 0, gg' > 0, C^ > 0) and (c = -1, J = 2, /3/3' > 0, AA' > 
0, gg' > 0, C^ > 0), respectively. Here, \h\ = 0.56, I5I = 0.75, l^'l = 0.25, |A| = 0.4546, |A'| = 0.2667, |/3| = 0.909, \P'\ = 0.533, 
Id = 0.727 and \va\ = 0.364. 
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FIG. 29: The bound state solution of the [D* D'i]°i system with (c = -1, J = 0, (3/3' < 0, AA' < 0, gg' < 0, (zu < 0), (c = 
-1, J = 0, ;g/3' < 0, AA' < 0, gg' < 0, (-w > 0), (c = -1, J = 0, /3/3' > 0, AA' < 0, gg' < 0, (zu > 0), (c = +1, J = 0, /?/?' < 
0, AA' < 0, gg' < 0, (uj < 0), (c = +1, J = 0, /3/3' < 0, AA' < 0, gg' < 0, (uj > 0), (c = +1, J = 0, 13/3' < 0, AA' < 0, gg' > 
0, Cra7 < 0), (c = +1, J = 0, /3/3' > 0, AA' < 0, gg' < 0, (w < 0) and (c = +1, J = 0, /3/3' > 0, AA' < 0, gg' < 0, (w > 0), 
which correspond to diagrams (a), (b), (c), (d), (e), (f), (g) and (h) respectively. Here, we take l/i. \h\ = 0.56, \g\ — 0.75, 
\g'\ = 0.25, |A| = 0.4546, |A'| = 0.2667, |/3| = 0.909, |/3'| = 0.533, |C| = 0.727 and \m\ = 0.364. 
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FIG. 30: The bound state solution of the [D* D'i]°i system with (c = -1, J = 1, /3/3' < 0, AA' < 0, gg' < 0, (u^ > 0), 
(c = +1, J = 1, /3/3' > 0, AA' < 0, gg' > 0, C^ < 0), (c = +1, J = 1, /3/3' > 0, AA' < 0, gg' < 0, (tx^ > 0), (c = +1, J = 
1, /3/3' > 0, AA' < 0, gg' < 0, C^^ < 0), (c = +1, J = 1, ^/3' < 0, AA' > 0, gg' < 0, Cw < 0) and (c = +1, J = 1, /?/?' < 0, AA' < 
0, gg' < 0, ("tJ7 < 0), which correspond to diagrams (a), (b), (c), (d), (e) and (f) respectively. Here, we use Ih. \h\ — 0.56, 
\g\^ 0.75, \g'\ = 0.25, jA] = 0.4546, |A'| = 0.2667, |/3| = 0.909, |/3'| = 0.533, |C| = 0.727 and \vj\ = 0.364. 
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FIG. 31: The bound state solution of the [D* D'i]°i system with (c = -1, J = 2, /3/3' < 0, AA' > 0, gg' > 0, (u^ > 0), 
(c = 4-1, J = 2, 13/3' > 0, AA' > 0, gg' > 0, C^ < 0), (c = +1, J = 2, /3/3' < 0, AA' > 0, gg' > 0, (tx^ > 0), (c = +1, J = 
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0, gg' > 0, ("tJ7 < 0), which correspond to diagrams (a), (b), (c), (d), (e) and (f) respectively. Here, we use Ih. \h\ — 0.56, 
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